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^ ' We show that homogeneous and isotropic cosmological models with radiation and scalar 

^ ' field have, in general, two possible chaotic exits to inflation. These models may describe 

^ the early stages of inflation. The central point of our analysis is based on the possible 

, existence of one or two saddle-centers points in the phase space of the models. 

^ . 

1 Introduction 

>' 

CO ' In the preceding paper we have examined the role of anisotropy in inducing chaotic 

. dynamics and chaotic exits to inflation for the Bianchi IX comsologial models. How- 

ever, chaotic behavior seems to be a general feature in the early satges of inflation. 
We ilustrate this in the realm of rather simple homogeneous and isotropic cos- 
' mological scenario, in which matter is represented by radiation and scalar fields. 

^\ . Non-integrability and chaos are invariantly characterized by the topology of homo- 

' clinic cylinderfB in the phase space. The relevant manifestation of chaos is that even 

small fluctuations of the scalar field are able to preclude or induce the universe to 



o 



escape to the inflationay phase with two distinct channels of exit to inflation. 
2 The Dynamics of the Model and the Chaotic Exit to Inflation 



The dynamics of closed homogeneous and isotropic cosmological models with scalar 
^ . field and radiation is governed by the following Hamiltonian 



H = -^^+^+ia-a^V{^)-^^0, (1) 



where a, Lp are the scale factor and the scalar field, respectively, whereas Pa, Ptp are 
the conjugated momenta. From the conservation of energy and the state equation 
p = p, it follows that p = Eq a~^, Eq being a constant of integration. V{(p) is the 
potential of the scalar field. 

There is one or more critical points in the finite region of the phase space 



denoted generically by P whose coordinates are gq — y jj-j, f — fo,Pa = Pip = 0, 

where (po is solution of the equation V'{(po) = and Ae/ = V{(po). The energy 
associated to the critical point is Eq ~ Ecr = foo- Assuming Agj > 0, this 
critical point represents the static Einstein universe. The stability analysis reveals 



the following eigenvalues associated to P: Ai_2 = ^ y | ^e/, A3^4 = ± ^y—V"{ipo). 
Two possibilities arise: if V"{(po) > 0, A3_4 are pure imaginary indicating that P 



1 



is a saddle-centerl, which induces a cylindrical topology in the phase space. These 
eigenvalues are associated to the hyperbolic and rotational motions around a small 
neighborhood of P. If l/"(<po) < 0, A34 are real with opposite signs meaning 
that P is a pure saddle. The phase space also has two critical points at infinity 
corresponding to the de Sitter solution, one acting as an attractor (stable de Sitter 
configuration) and the other as a repeller (unstable de Sitter configuration). The 

scale factor a{t) approaches to the de Sitter attractor as a ~ ev^^^*, whereas 
Pip — Q,^p — Lpi). Finally, another very important feature of the resulting dynamical 
system is the admittance of invariant manifolds AA defined generically hy = 
fOj Pip = 0, where the the dynamics is governed by a two-dimensional system. Such 
a system is integrable and its integral curves represent homogeneous and isotropic 
universes with radiation and cosmological constant A^f (cf. Fig. 1). 




Figure 1: Integral curves on the invariant man- Figure 2: Chaotic exit to inflation: 30 orbits 
ifold = 0, ip = ipo for 7 = 4/3 (radiation). taken in a sphere of radius R = 10~^ about a 

point ,3(1 on the separatrix. 

The numerical experiments are performed using the package Poincar^. Two 
distinct situations are considered: (a) V{lp) = A + ^m'^ip'^, and (6) V{(p) = A + 
J (</j^ — cr^)^, where the cosmological constant A is included, m is the mass of 
the scalar field, A and a are positive constants. For the first choice, the phase 
space has only one critical point P identified as a saddle-center. There is also 
one invariant manifold A4. In the same way, as pointed in Ref. [3], the phase 
space under consideration is not compact, and we will actually identify the chaotic 
behavior associated to the possible asymptotic outcomes of the orbits in this phase 
space, namely, escape to de Sitter state attractor at infinity (inflationary regime) 
or collapse after a burst of initial expansion. 

Let us consider initially case (a) . The behavior of the orbits near the separatri- 
ces that belong to the invariant manifold Ai is performed numerically. We remark, 
however, that the separatrices define the regions of collapse and expansion in M, 
but not in the full 4-dimensional phase space. Let us consider a point belonging 
to the separatrix. Around this point, we construct a 4-dimensional sphere of initial 
conditions in the phase space with arbitrary small radius R, say R = 10-4, that 
represent expanding models after the initial singularity. The main result is showed 
in Fig. 2, where there is always a domain of energy such that part of the orbits es- 
cape to the de Sitter phase and another part collapse, resulting in an indeterminate 
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outcome. Note that, before the escape or collapse, the orbits visit a small neigh- 
borhood of the saddle-center P from which the petition of the total energy into 
rotational and hyperbolic modes about P is chaotifJ. Putting in another words we 
state that the boundaries of collapse and escape are chaotically mixed. 

The choice of the potential (b) with symmetry breaking offers a richer phase 
space dynamics. Three critical points, two of them are saddle-centers whereas the 
remaining is a pure saddle, coexist. Also, three invariant manifolds, each containing 
one critical point are present. Repeating the numerical experiments of taking a 
small sphere of initial conditions very close to the separatrices (or any other curve) 
of one of the invariant manifolds for which the critical point is a saddle-center, we 
observe the same chaotic behavior with respect to collapse and escape to the de 
Sitter phase, as indicated in Fig. 2. However, if inital conditions are taken near the 
invariant manifold of the pure saddle, we note that, depending on the total energy 
Eq, two chaotic escapes to inflation are possible. In the first, there is always domain 
of energy Eq such that all orbits evolve towards one of the saddle-centers, perform 
some oscillations around it to proceed to collapse or escape to the de Sitter phase as 
showed in Fig. 3. In the second, all orbits, after visiting one of the saddle-centers, 
return to an infinitesimal neighborhood of the pure saddle, where part of them 
collapse and another escape to inflation (see Fig. 4). In this case, the partition of 
the total energy into the hyperbolic motions around the pure saddle is completely 
chaotic. 




Figure 3: Chaotic exit to inflation after oscilla- Figure 4: Another chaotic exit to inflation: the 
tions about the saddle-center. orbits return to the pure saddle before collapse 

or escape. 
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